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Abstract 

Classical nucleation theory is used to estimate the free-energy barrier to nucleation of the solid phase of particles 
interacting via a potential which has a short-ranged attraction. Due to the high interfacial tension between the fluid 
and solid phases, this barrier is very large, much larger than in hard spheres, ft is divergent in the limit that the 
range of the attraction tends to zero. We predict an upper limit on nucleation in good agreement with the results of 
experiments on the crystallisation of proteins. 



1 Introduction 

It is difficult to crystallise globular proteins yet it is impor- 
tant to do so as X-ray crystallography can be performed 
on crystalline proteins and thcLfiill structure of the protein 
in its native state determined.Era Knowing the structure of 
a protein is crucial to understanding its function. Proteins 
are of course vital to all life. Crystallising a protein is diffi- 
cult as there seem to be only a few conditions under which 
a protein solution crystallises but many under which the so- 
lution is cither stable on all accessible time scales or useless 
amorphous aggregates form. This is in sharp contrast to 
the simpler system of hard-sphere-like colloids. The crys- 
tallisation of these is routine. Here we will examine the 
first stage of crystallisation: nucleation. Wc will do so us- 
ing classical nucleation theory. Although not a quantitative 
theory, it has the merit of simplicity and incorporates the 
essential features of nucleation. These are that nucleation 
is a large amplitude, localised, fluctuation which forms a 
microscopic crystallite of the crystalline phase. 

Because of their phase diagrami'i it has been sug- 
gestedu that the interaction between globular protein 
molecules can be approximated by a steeply repulsive core, 
which prevents two protein molecules interpenetrating, and 
an attraction with a range that is much less than the diame- 
ter of the protein; the surfaces of two protein molecules have 
to be within a tenth or less of their diameter for there to be 
a significant attraction. Here we will study a simple model 
potential of this type. In fact, the interactions between 
protein molecules are complex and include quite specific, 
non-spherically symmetric interactions,!^ thus it should be 
borne in mind that our model potential is extremely sim- 



plified. 

Classical nucleation theory requires as input the chem- 
ical potential of the metastable fluid phase from which 
the crystal forms, the chemical potential of the crystalline 
phase and the surface tension. The last of these three, the 
surface tension, often poses a problem as the surface ten- 
sion is known for only a very few model potentials. Indeed 
even for the simplest potential, the hard-sphere potential, 
accurate values foiithe surface tension have only been cal- 
culated recently .Lm However, recently the author has cal- 
culated the surface tension for a potential consisting of a 
hard core and a very short-ranged attraction.El 

In the next section we will define our model potential 
and describe its phase behaviour. We will give the required 
expressions for the chemical potential in the fluid and solid 
phases. Then in the third section we will obtain the classi- 
cal nucleation theory predictions for our model. The fourth 
section compares our theory with experiment and the fifth 
is a conclusion. 



2 Model and bulk phase diagram 

First, we define the hard-sphere+square-well potential]! It 
is the spherically symmetric pair potential u(r) defined by 
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where a is the hard-sphere diameter, and r is the separation 
between the centres of the spheres. The potential depends 
on one parameter, 5, the ratio of the range of the attrac- 
tion to the hard-sphere diameter. We will always consider 
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this range to be small: 6 -C 1. The first person to consider 
very short-ranged attractions was Baxtem who considered 
a potential with zero range, 5 = 0, and with a well-depth 
e/kT adjusted so that the second viriai coefficients was of 
order unity. T is the temperature and k is Boltzmann's 
constant. This model is often termed the sticky-sphere 
model. Within it the second viriai coefficient is used as 
a temperature-like variable. However, Stel£3 showed that 
the sticky-sphere model was pathological, its fluid phase is 
unstable at all non-zero densities. Therefore, we will not 
consider this model but instead will follow Stcll pziipa we 
take the limit 5 — ► 0, thus obtaining his (3q modelErllj 

The bulk phase behaviour of the j3o model is described 
in Rcfs. H and|lj, and is displayed in Fig. |. If 

the double limit 6^0 and kT jt — > is taken such that 
T > T co u then the equilibrium phase behaviour is identical 
to that of hard spheres. If T < T co u then the behaviour 
is radically different: a close-packed solid coexists .with, 
fluid phase of zero density. The temperature T co u isE3'Or 



kT coU 
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The close packed density is 2 1 / 2 cr~ 3 which corresponds 
to a volume fraction r\ occupied by the hard cores of 
rj = 2 1 / 2 7r/6 ~ 0.74. At close packing there is an isostruc- 
tural solid-solid transition.113 The bulk phase behaviour 
of hard spheres with a short-range square-well attraction, 
0.06 > 6 > 0, is qualitatively similarly It is essentially 
Fig. ^ but with the corners rounded off and the horizontal 
(vertical) lines no longer horizontal (vertical) but sloping, 
see Refs. |l3 , 15 and [l6|. The coexisting densities no longer 
change discontinuously but only rapidly near T co u. For 
S > 0.06 there is no solid-solid transition but there is no 
vapour-liquid transition until considerably larger values of 
8. 

2.1 Phase coexistence for nonzero 5 

At temperatures significantly below T co u, the density of the 
fluid phase which coexists with the solid phase is very low 
(zero in the 8 — > limit). So, we will approximate the 
chemical potential of the fluid phase by that of an ideal 
gas. Around the temperature T co u the density of the fluid 
phase which coexists with the solid rapidly increases (in the 
5^0 limit it is horizontal, see Fig. ||). So the accuracy of 
the approximation will rapidly worsen as the temperature 
reaches T co u . The chemical potential, pf, of an ideal gas 



fj,f = fcTlnp, 



(3) 



where p is the fluid's number density of spheres times 
er 3 to render it dimensionless. We have neglected a term 
ln(A/cr 3 ), where A is the integral over the momentum de- 
grees of freedom. 

At low temperature, the Helmholtz free energy per par- 
ticle, a s , of the solid phase can be calculated by using a 



cell theory ,0 and assuming that the attractions are strong 
enough so that they force the lattice constant to be small 
enough so that each sphere is within the range of the atj 
tractive interaction of all twelve of its nearest neighbours £3 
This fixes the energy per particle at — 6e and the lattice 
constant at <r(l + cS), where c is a constant a little less 
than one, sej-Rdj. |l4] . Then the Helmholtz free energy per 



particle 



(4) 



where qp is the partition function of a single particle 
trapped in a cage formed by its twelve neighbours fixed 
at their points in a perfect lattice. This lattice is either 
face-centred cubic or hexagonal close-packed, in each a par- 
ticle has the same arrangement of nearest neighbours so it 
does not matter which. As the lattice constant is suffi- 
ciently small, throughout the cage the particle attracts all 
twelve neighbours and so its energy is constant at — 6e. 
The partition function is then just the volume available to 
the centre of mass of the particle times A -1 exp(6e//cT), 
where A -1 is the integral over the momentum degrees of 
freedom. The volume is of the order of (cSa) 3 , i.e., the 
cube of the difference between the lattice constant and the 
hard-sphere diameter. The particle at its lattice position 
is <r(l + cS) away from each of its neighbours and cannot 
move within a of them due to the hard-sphere interaction. 
Thus, it can move a distance cSa in the direction of each 
of its neighbours. Actually, the volume to which the par- 
ticle is restricted is quite complicated, see Ref. but ^ 
is not far from (c5a) 3 . We will approximate it by (Scr) 3 ; 
an approximation which becomes increasingly accurate as 
S decreases £-3 Then, we have for qp, 



q P = (<5o-) 3 A- 1 exp(6e/fcT) 
Inserting this expression for qp into Eq. (0 
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where we have neglected a term ln(A/cr 3 ). The chemical 
potential p s = a s + p s / p where p s is the pressure and p is 
the density. At coexistence the pressure is very low, note 
that the gas density is zero in Fig. [j]. Therefore, p s /p con- 
tributes a negligible amount to the chemical potential for 
the /3o model. This enables us to equate a s and fi s as we 
have done in Eq. (g). For the /3 model all temperatures 
below T co u in the d —> limit are low in the sense that 
the solid is close-packed and so no different from the solid 
at zero temperature. When 6 is nonzero the accuracy of 
our low temperature approximation, Eq. (g), will rapidly 
worsen as T co u is approached from below. 

Equating the chemical potentials of the fluid, Eq. (||), 
and solid, Eq. (^|), phases yields an equation relating the 
density of the coexisting fluid to the temperature of coexis- 
tence, at a given value for 6. We can rearrange the equation 
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to give the temperature of coexistence, T co , as a function 
of the density of the coexisting fluid phase, pf, 
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each value of B2 or r to a temperature T and vice-versa. 
However, if 8 = then B2 is not a well-behaved function 
of temperature. If we take the double limit 8, kT/e — > 
thenfca 

B2 _ ( 1 T > Tb20 /, o\ 



As we cool a fluid phase of density pf, then, at equilibrium, 

the solid phase appears at the temperature T co of Eq. (Q). where Tb2o is given by 
Using Eq. (0) for T co ;/ we can express T co as a fraction of 
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T coll l-ln P/ /[31n(V*)]" 

As 8 is small the ratio T co /T co u is near one for a very wide 
range of p$. For example, when <5 = 0.05, T co /T coll — 0.80 
and 0.66 for pj — 0.1 and 0.01, respectively. 

2.2 The second virial coefficient 

For the square- well potential of Eq. ([!]) with 8 <C 1, the 
second virial coefficient is 



B 2 = B' 2 ls - 2 7 rer i (5(exp(e/fcT) - 1), 



(9) 



where B' 2 ls = (2ir/3)a 3 is the second virial coefficient of 
hard spheres. Generally we will be working at sufficiently 
large values of e/kT that we can neglect the —1 term in 
parentheses in Eq. (^). Therefore, we will neglect this 
term from now on. From the second virial coefficient Bax- 
ter defined the t parameters! 



(10) 
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which has been described as an effective temperature al- 
though it is not a simple reduced temperature as it is not 
a simple ratio of the thermal energy, kT, to an energy. It 
can be written directly in terms of the temperature and the 
range 8 

T = 12(5exp(e/fcT)' ^ 

Straightaway, we can find r at the temperature T co u, t co u, 
by substituting Eq. (|) into Eq. (|ll|), 
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thus as 8 tends to zero, t co u diverges. t co u is the value 
of t at Tcou, when the fluid-solid transition widens from 
the relatively narrow coexistence region of hard spheres to 
coexistence between a dilute gas and a near close-packed 
solid. 

The experiments measure the strength of the attractive 
interaction by measuring the second virial coefficient via 
scattering whereas we measure it via the ratio of the well- 
depth to the thermal energy, e/kT. Thus we consider here 
the problem of relating B 2 and r to e/kT. For nonzero 
8, B2 and hence x is a well-behaved single- valued func- 
tion of temperature and so it is possible to uniquely map 



ln(l/«J) 
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All temperatures above Tb20 map onto a second virial co- 
efficient equal to B2 8 , and all values below Tb 2 o ma P onto 
—00. Conversely, all nonzero and non- infinite values of 
(B2 — B2 S )/B 2 S , and hence of r, correspond to the tem- 
perature Tb2o, see the discussion of Stell.t2l Which is why 
t is useless to describe the phase diagram of spheres with 
an attraction of zero range. 

Away from the double limit, where 8 is small but 
nonzero, the variation of r with temperature depends 
strongly on the range 8. We plot r, Eq. (pi]), against the 
reduced temperature kT/e in Fig. 0. At any given reduced 
temperature, kT/e, t varies as 8 . Determining e from 
scattering data for r thus requires an accurate knowledge 
of 8. This is not available for globular proteins, see section 
4. 

2.3 The metastable fluid phase 

Nucleation occurs from a metastable phase, i.e., a phase 
which is not the equilibrium phase. We cannot obtain 
the nucleation rate in a metastable phase without mak- 
ing assumptions about this phase. We will assume that 
the metastable phase is close to an ideal gas. Justifica- 
tion for this assumption is provided by the fact that the 
densities we will consider will be quite low p < 0.1 and 
the second virial coefficient, Eq. (^), will not be too large, 
B2 1 B^ of order unity. The second virial coefficient is the 
coefficient of the leading order correction to the ideal gas 
pressure in a density expansion for the pressure J3 How- 
ever, assuming ideal gas behaviour will not be correct if 
the attractive interactions are strong enough to cause non- 
equilibrium effects such as gelation or in the vicinity, of a 
metastable fluid-fluid transition 1*3 .Both gelationEHEl and 
a metastable fluid-fluid transitionEJ have been observed in 
experiment. 



3 Classical nucleation theory 

Classical nucleation theory&@ assumes that the nucleus of 
the solid phase has a free energy AF which is the sum oLa 
bulk term and surface term. See the book of DebenedettiE3 
for an excellent introduction to classical nucleation the- 
ory. The bulk term is equal to the number of spheres 
in the nucleus, n, times the chemical potential difference 
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Ap = fi eq — jj, m , where p eq is the chemical potential of the 
coexisting equilibrium fluid and solid phases, and fi m is 
the chemical potential of the metastable fluid phase which 
contains the nucleus. Our Eq. (^|) gives the chemical po- 
tential of the solid phase at coexistence p s which is equal 
to Heq . The surface term is the surface area of the nucleus 
times the surface tension 7 of the bulk interface between 
the coexisting solid and fluid phases, i.e., any difference 
between the spherical interface between the solid nucleus 
and the surrounding metastable fluid phase, and the flat 
interface between the coexisting solid and fluid phases is 
ignored. The interfacial tension of an interface between a 
solid phase and a fluid phase depends on the orientation of 
the interface with respect to the lattice of the solid. To get 
around this problem, either some orientationally averaged 
interfacial tension or the interfacial tension with the inter- 
face in the plane of a low index lattice plane are used for 
7. The surface area of the nucleus is obtained by assum- 
ing that the nucleus is a perfect sphere of radius r which 
is related to n by assuming that the density of spheres 
within the nucleus is equal to the bulk density of the solid 
phase. For the 0q model this is the close-packed density .of 
hard spheres 2 1 / 2 ct~ 3 . When 5 is nonzero the density isli-3 
2 1 / 2 Il + c<5)~ 3 er~ 3 , where c is a constant which is less than 
onc.li-3 We will see below that its value if irrelevant and so 
we do not attempt to estimate it. Then the surface area of 
the nucleus is n 2 / 3 (187r) 1 / 3 (l + cS) 2 a 2 . So, 

AF(n) = nAp + n 2 / 3 (187r) 1 / 3 7 (l + c8fa 2 (15) 

The first term in AF is the bulk term, which is negative 
and decreases linearly with n, and the second term is the 
surface term which is positive and increases as n 2 / 3 . Thus 
AF passes through a maximum, denoted by AF* , at some 
value of n, denoted by n*. The nucleus with n* spheres 
is called the critical nucleus. This is the nucleus of the 
solid phase at the size at which its free energy is greatest. 
As the frequency with which a fluctuation occurs scales as 
exp(— free energy cost/fcT), the formation of the fluctua- 
tion with the highest free energy cost is the rate limiting 
process. Thus, classical nucleation theory predicts that the 
nucleation rate is proportional to exp(— AF*/kT) where 
AF* = AF(n*) is the maximum in the free energy of Eq. 
(H). Taking the derivative of Eq. (||) with respect to n 
and equating to zero yields the number of spheres in the 
cluster at the top of the barrier 
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and inserting this value of n in Eq. (jig) 



AF* 



8^ (7(1 + c<5) 2 f7 2 ) : 
3A/I 2 



(16) 



(17) 



Now, the interfacial tension, 7m, for an interface in 
the 111 crystal plane, between the coexisting close-packed 



solid phasa,and infinitely dilute gas phase of the /3q model 
is given byu 

7lll = 3 1/2 ecr- 2 6 = 0. (18) 

Which, like our expression for the chemical potential, Eq. 
@, is a low temperature approximation also valid when 
< 5 <C 1, but which breaks down, when 5 is nonzero, 
as T co u is approached from below. Other lattice planes 
have slightly higher tensions, e.g., the 100 plane has a sur- 
face tension 7100 = 2ea~ 2 . We will approximate 7 for all 
small values of 6 by the surface tension of the 111 plane 
7m, which will vary with 5. The tension 7m is the en- 
ergy penalty due to the broken bonds at the surface of the 
solidE Therefore, it scales as one over the square of the 
lattice constant, (1 + c8)a. So, for all small values of 5 we 
have 



jjl + cSfa 2 _ 1/2 _ 
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kT 



kT 



ln(l/5) 



T 



(19) 



where the second equality was obtained using Eq. In 
Rcf. fj] it was shown that the gas-solid interface is narrow, 
only one sphere wide, and so the assumption_pf_classical 
nucleation theory that the interface is narrowErEil is fully 
justified. 

Then, the free-energy barrier to nucleation, AF*, ex- 
pressed as a function of temperature and the excess chem- 
ical potential over that at equilibrium, is 



AF* _ 3 1 / 2 7 r[ln(l/ ( 5)] 3 (T/T coll ) 
kT (Ap/kT) 2 

The corresponding expression for n* is 

, _ 3 1 / 2 2^ [ln(l/<5)] 3 (T/T coll ) 
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Equations © and © are valid below and not too close 
to Tcou. From them we see that at fixed chemical potential 
difference, Ap/kT, and ratio of the temperature to T co u, 
as 6 decreases both the free-energy barrier and the size of 
the critical nucleus increase as [ln(l/(5)] 3 . The increasing 
surface tension, Eq. (|l^), as S decreases is increasing the 
barrier to nucleation. 

We would like to express AF* and n* purely in terms 
of temperature, density of the metastable fluid phase, pf, 
and 6. To do this we need to know A/i = fj, s — fi m . Now, 
[j, a is given by Eq. (||). This leaves the effective chemical 
potential of the metastable fluid phase. We will estimate 
this by assuming that the metastable fluid phase closely 
resembles an ideal gas, see Section 2.3. This is an accurate 
approximation near coexistence. So, with this assumption, 
fj, m is given by Eq. (||) with p = pf, and so A/i is 



(22) 
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If we use Eq. (Q) to obtain a value for T co at the given value 
for pf then we can define a reduced temperature t = T/T co ; 
a function of T and pj. The reduced temperature t is the 
temperature as a fraction of the coexistence temperature 
at the same value of the fluid density. In terms of t, Eq. 
(|22]) can be rewritten as 



Ap 
"fcT 



31n(l/<5)[l - r 
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31n(l/<J) 



(23) 



This is the excess chemical potential of a metastable fluid 
phase at a density pf and at a fraction t of the tempera- 
ture at which a fluid of this same density pf coexists with 
the solid phase. Ap is of order kT when (T co — T)/T co = 
1 — t ~ 1/ ln(l/<5). To achieve a chemical potential in the 
metastable fluid which is of order kT higher than at equi- 
librium we need only go a small fraction of T co below T co . 
When the temperature difference (T co — T)/T co is finite, 
Ap diverges in the i5 — ► limit. 

We can substitute Eq. (23) for Ap, and Eq. (||) for 
Tcoii into Eq. (||^) to obtain the free energy barrier to nu- 
cleation as a function of temperature and the range of the 
attraction, <5, 



AF* _ 3- 3 / 2 7rln(l/<5) (1 - ln P/ /[3 Hl/5)]) 
kT ~ *(1 — i) 2 ' 



(24) 



Which is the free energy barrier to the nucleation of the 
solid phase in a metastable fluid phase at a density p / and 
at a fraction t of the temperature at which a fluid of this 
same density pf coexists with the solid phase. At a given 
fraction, t, of the phase transition temperature, the free en- 
ergy barrier to nucleation increases as ln(l/5), as the range, 
<5, decreases. The corresponding expression for n* is very 
simple, 

3-5/2 27r 

"* - TTW- (25 » 

which is not a function of the range 5 at fixed t. It is the 
number of spheres in the critical nucleus in a metastable 
fluid phase at a density pf and at a fraction t of the temper- 
ature at which a fluid of this same density p / coexists with 
the solid phase. Our Eqs. ( p4| ) and ( p5| ) for AF* and n*, 
respectively, rely on our approximations for 7 and Ap be- 
ing accurate, which they are provided we are not too close 
to T co u when 6 is nonzero, and provided the metastable 
phase is not too unlike an ideal gas. Also, the free energy 
is derived by assuming that we can treat it as the sum of 
a bulk part and an interfacial part. This is reasonable if n 
is large but breaks down when it is small. Thus, classical 
nucleation theories break down when n* drops below about 
twenty. From Eq. (||), n* = 20 when t = 0.73. Thus be- 
low t ~ 0.7 the critical nucleus is so small that classical 
nucleation theory can no longer be expected to be reliable. 

The free energy barrier AF* /kT and the number of 
spheres n* at the top of the barrier are plotted in Figs. 
H and ||, respectively, for three different ranges 6 = 0.1, 



0.05 and 0.01. At any value of t the shorter the range 
the larger is AF* /kT. This is true even though as the 
range decreases Ap/kT, Eq. (p3|), increases. AF* /kT 
varies as (7<r 2 /kT) 3 / (Ap/kT) 2 and as both the surface ten- 
sion 7 and Ap/kT vary as ln(l/<5), then AF* /kT varies as 
ln(l/£), Eq. (|24|). But the number of spheres at the top 
of the barrier is independent of range, see Fig. ^. Thus we 
predict that when the attraction is very short-ranged, the 
free energy barrier to nucleation is larger than it is when 
the attraction is longer-ranged but the number of spheres at 
the top of the free energy barrier, n* , is the same. As both 
the free-energy barrier and Ap/kT increase with decreas- 
ing range then the free-energy barrier will only be small 
enough for nucleation to occur when Ap/kT is very large. 



4 Comparison with protein nucle- 
ation 

We have specified the state of our model system by the den- 
sity p and the reduced temperature kT/e, which is the ratio 
of the thermal energy to the attractive energy. For protein 
solutions the attraction is a free energy not an energy, and 
so its strength e is a function of temperature, salt concen- 
tration and other variables. The dependence on tempera- 
ture is then rather complex. For example adding polymer 
induces a depletion attractioncj'Ej which is proportional to 
temperature and to the concentration of polymer. All this 
makes a direct determination of the ratio of the thermal 
energy to the attraction free energy all but impossible in 
experiment. What can often be measured in experiment is 
the second viral coefficient. George and WilsonEa use the 
second virial coefficient B2 as a measure of the strength of 
the attraction, and Rosenbaum et alB use Baxter's r pa- 
rameter. The second virial coefficient, Eq. (^), is a function 
of e/kT and 6. Thus to convert a temperature into r or vice 
versa we need to know the range of the attraction between 
protein molecules. We do not know this range, also we 
do not know the extent to. which it varies from protein to 
protein. Rosenbaum et al.u study lysozymc, effective diam- 
eter 3.4nm and STA, effective diameter 1.3nm. Only if the 
range of the attraction scales with the size of the protein 
will the range of the attraction be the same in each case. In 
the absence of an estimate for the range of the attraction 
between protein molecules we have performed calculations 
for a few values of the range, S. 

Rosenbaum et aiH plot the values of Baxter's r pa- 
rameter for which a number of globular proteins crys- 
tallise. They conclude that "crystallising conditions occur 
for 0.06 < t < 0.15". In order to compare the results of 
classical nucleation theory with experiment we rewrite the 
free-energy barrier, Eq. (0), in terms of r. To do this we 
use Eq. @ to substitute r for e/kT in Eqs. © and (||) 
and then we substitute the resulting expressions for 7 and 
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A/z into Eq. ( Jl7| ) to obtain 

AF* _ 3 1 / 2 8 7 r|lnl25r| 3 
kT ~ (3 ln(l/<5) - In p f + 6 In 12£r) 2 ' ^ ' 

We have plotted AF* /kT as a function of r in Fig. ||. 

For S = 0.1 and 0.05, the free-energy barrier is a very 
sharply varying function of t. A little below t = 0.15 
and 0.2, for 6 = 0.1 and 0.05, respectively, the free-energy 
barrier increases to huge values. The sharply varying free 
energies and the fact that the nucleation rate depends ex- 
ponentially on the free-energy means that there will be a 
well defined upper limit to the value of r at which nucle- 
ation occurs during a given time scale. Our prediction, for 
very reasonable ranges d — 0.05 and 0.1, of an upper limit 
to nucleation at r around 0.15 is in good agreement with 
experiment. Despite the crudity of our theory this agree- 
ment may not be fortuitous. We do not know the range 
S that best corresponds to the interaction potential of any 
globular protein and our calculation of the free energy bar- 
rier may be out by tens of kT. However, when the range 
is of the order of a few % of the hard-sphere diameter then 
changing the range by a factor of two or more and adding or 
subtracting several tens of kT to AF* do not alter our pre- 
diction that nucleation cannot occur above r = 0.15 ±0.05. 
Because of the nucleation rate's exponential dependence on 
the distance into the metastable regime it is common to ob- 
serve a well-defined threshold beyond which nucleation is 
not observed. BEj Thus our finding of a well-defined upper 
limit on nucleation is no surprise. However, the fact that 
this upper limit occurs at a value of r which depends very 
weakly on the range of the attraction, 5, is more interesting. 
It is a possible explanation of why many different proteins 
have approximately the same upper limit even though the 
effective interaction potential may differ from protein to 
protein. 

If we compare Figs. || and || we see that when the bar- 
rier height is plotted as a function of t, the shorter the 
range the steeper the curve followed by the barrier height, 
whereas when the barrier height is plotted as a function of 
t it is the other way around. Changing variables from t to 
r reverses the trend in the change of the free-energy bar- 
rier with S. This is a consequence of the increasingly rapid 
variation of r with temperature as 5 decreases. It suggests 
caution in interpreting r as a temperature-like variable. 

Although the interaction between globular proteins is 
complex and poorly understood, thprp is an experimental 
system with a known, short range.EBEj This is hard-sphere- 
like colloidal particles with added non-adsorbing polymer 
with a radius of gyration much less than the diameter of 
the particles. Then the polymer induces a depletion at- 
traction!^} E3 with a range S approximately equal to the 
radius of gyration diidded by the diameter of the colloidal 
particles. Ilett et a/.E£lEj studied a system with a range S 
of approximately 0.08. However, as the coexistence den- 
sities were not determined accurately we cannot compare 



our theoretical predictions with these experiments. Qual- 
itatively, they find that crystallisation only occurs if they 
are not too far into the metastable regime, if t is not too 
small. If t is too small then a gel forms. As for proteins 
there is a window within which crystallisation takes place. 

5 Conclusion 

We have calculated the free-energy barrier to nucleation, 
within classical nucleation theory, for spherically symmet- 
ric potentials consisting of a hard-core and a short-range 
attraction. A simple model potential for globular proteins 
and some colloidal systems. The short-range attraction 
makes nucleation more difficult: the surface tension, Eq. 
dl9|), is larger than for many other systems, e.g., hard 
spheres, which means that at the same value of Afi/kT 
the free energy barrier to nucleation is higher when there 
is a short-ranged attraction than for hard spheres. The 
surface tension for-hard spheres is, from density- functional 
theory, 0.3fcTcr _2 ,EHa whilst with an attraction of range 10% 
of the diameter of the particle, <5 = 0.1, at a temperature 
(2/3)T co h, it is 5.9kT<r~ 2 . The surface tension for hard 
spheres is an order of magnitude smaller and the nucle- 
ation rate is very sensitive to the surface tension, it varies 
as exp [-(7cr 2 /fcT) 3 ] at constant Afi/kT. 

The large barriers to nucleation predicted by classi- 
cal nucleation theory are certainly consistent- with the 
difficulty experienced in crystallising proteinsmSEa and 
the non-equilibrium behaviour observed ibr colloidal par- 
ticles with a short-ranged attraction£3'c2l However, the 
uncertainty surrounding the interaction between protein 
molecules makes comparison with theory, which is for a spe- 
cific model potential, very difficult. Expjemnental results 
for the location of the onset of nucleationEBEJ for colloidal 
systems for which the the interaction potential potential 
and phase coexistence boundaries are both known would 
be very useful in understanding protein crystallisation and 
in testing theoretical predictions. 

The barrier to nucleation decreases as the temperature 
drops more and more below the transition temperature T co , 
i.e., as t decreases, and we move further into the metastable 
region of the phase diagram. Thus we can lower the free en- 
ergy barrier by moving deeper into the metastable regime. 
However, by doing so we may encounter difficulties such 
as the. ff-wd. becoming sufficiently cold to form a gel-like 
stateEj'E3Ej in which the dynamics are at least partly ar- 
rested, preventing the ibrmation of the crystalline phase. 
Ten Wolde and Frenkellia have suggested that a metastable 
fluid-fluid transition may offer a way round this problem. 

An interesting prediction of classical nucleation theory 
is that the critical nucleus at a given value of the free-energy 
barrier shrinks as the range of the attraction, 5, decreases. 
Thus under the conditions in which the free-energy bar- 
rier is small enough to permit nucleation on a reasonable 
time scale the nucleus becomes smaller and smaller as 5 
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decreases. Indeed it becomes so small that classical nucle- 
ation theory breaks down, even when AF* is still large, of 
order tens of kT . 
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Figure 1: The phase diagram of the /3q model.ll^rBEJ The solid lines denote the coexisting densities at the fluid-solid 
and solid-solid transitions. The density difference between the coexisting phases at the solid-solid transition is zero on 
this scale.113 The letters denote the fluid (F) and solid (S) phases, and the two-phase region of fluid-solid coexistence is 
marked with the number 2. The dashed lines are tie lines connecting coexisting fluid and solid phases. 
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Figure 2: Baxter's r parameter, Eq. (JlT|) , plotted as a function of the ratio of the thermal energy to the well depth, 
kT/e. The solid, dashed, and dotted curves are for ranges 6 = 0.1, 0.05 and 0.01, respectively. kTs^ajt = 0.43 and 0.22 
for S — 0.1 and 0.01, respectively. 

2.5 



1.5 



0.5 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.! 

kT/e 



0.9 1 



Figure 3: The free energy barrier AF*/kT is plotted as a function of t = T/T co for three ranges: 5 = 0.1 (solid 
curve), 5 = 0.05 (dashed curve) and 5 = 0.01 (dotted curve). For S = 0.1, 0.05 and 0.01, kT co /e = 0.65, 0.53 and 0.37, 
respectively. 
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Figure 4: The number of spheres at the top of barrier, n* , as a function of t = T/T co . The one curve plotted is for any 
(short) range S. 
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Figure 5: The free energy barrier AF* /kT is plotted as a function of Baxter's r parameter for three ranges: 5 = 0.1 
(solid curve), 6 — 0.05 (dashed curve) and S = 0.01 (dotted curve). All three curves are plotted down to r = 0.1 but in 
the case of S = 0.01 the size of critical nucleus, n* , drops below twenty at r = 0.2 and so the dotted curve from r = 0.1 
to 0.2 is really out of the range of applicability of classical nucleation theory. 
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